Formation of quarkonium states at RHIC 

R L Thews, M Schroedter and J Rafelski 

Department of Physics, University of Arizona, Tucson, AZ 85721, USA 

o 

C ~ J ■ Abstract. At RHIC the cross section for ce production will be large enough 

f^*^ ■ such that approximately 10 pairs will be produced in each central collision. If a 

^vj ' region of deconfined quarks and gluons is subsequently formed, one would expect 

that the mobility of the charm quarks will enable them to form J/ip through "off- 
t-*H, diagonal" combinations, involving a quark and an antiquark which were originally 

produced in separate incoherent interactions. We present model estimates of this 
<yj • effect, which indicate that the signal for deconfinement at RHIC may possibly be 

1^^ J/ip enhancement rather than suppression. 

> 

(O 1. Introduction 

A decrease in the number of observed J/ip in heavy ion collisions due to the screening 
of the color confining potential was proposed many years ago M as a signature of 
a deconfined phase. It is argued that as the system cools and the deconfined phase 
disappears, these heavy quarks will most likely form a final hadronic state with one 
of the much more numerous light quarks. The result will be a decreased population 
of J/ip relative to those formed initially in the heavy ion collision. 

Here we study a scenario which can only be realized at RHIC (and LHC) energies, 
Jj~H- where the average number of initially-produced heavy quark pairs Nq is substantially 

above unity in each central collision. Then one can amplify the probability of J/ip 
formation by a factor which is proportional to Nq , if and only if a space-time region 
of deconfined quarks and gluons is present. Realization of this result will depend on 
the efficiency of this new formation mechanism during the deconfinement period. We 
have developed a simple model to estimate the magnitude of this effect, and examined 
the sensitivity of results to various input parameters and assumptions. 

2. Suppression Factor 

For expected conditions at RHIC, almost all of the directly-produced J/ip will be 
dissociated even in peripheral collisions. To include the effects of our new formation 
mechanism, we parameterize the final J/ip number in each event as follows: 

Of the Nq charm quark pairs initially produced in a central heavy ion collision, let 
N\ be the number of those pairs which form J /ip states in the normal confining vacuum 
potential. At hadronization, the final number Nj/^, will contain a small fraction e of 
the initial number N\. The majority of Nj/^ will be formed by this new mechanism 
which we expect to be quadratic in the remaining Ao — N± heavy quark pairs, with a 
proportionality parameter (3. (We include in the new mechanism both formation and 
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suppression effects, since they occur simultaneously in the deconfined region.) The 
final population is then 

N J/i! =eN l +f3{No-N 1 ) 2 . (1) 

For each No initially-produced heavy quark pairs, we then average over the 
distribution of N±, introducing the probability x that a given heavy quark pair was 
in a bound state before the deconfined phase was formed. (This factor includes the 
effect of interactions with target and projectile nucleons). We finally average over the 
distribution of No, using a Poisson distribution with average value No, to obtain the 
expected < Nj/^, > final population per collision, 

< N JH> >= xN a (e + p(l - x)) + No(N + 1)0(1 - xf . (2) 

The bound state "suppression" factor Sji$ is just the ratio of this average 
population to the average initially-produced bound state population per collision, 

xN Q . 

Sj^=e + p(l-x)+l3 {1 ~ x) (JVb + 1) (3) 

x 

Without the new production mechanism, (3 — and the suppression factor is 

Sj/ip = e < 1. (Even the fitted parameter e contains some effects of the new 

mechanism, since formation can reoccur subsequent to the dissociation of an initial 

J/tp. Here we use it as an upper limit with which to compare the complete result.) 

However, it is possible that for sufficiently large values of j3 and No this factor 

could actually exceed unity, i.e. one would predict an enhancement in the heavy 

quarkonium production rates to be the signature of deconfinement! We thus proceed 

to estimate expected (3- values for J/ip production at RHIC. 

3. Model for J/tp Formation 

This model is adapted from our previous work on the formation of B c mesons B. 
Initial results for the J/ip application are found in Reference Q. For simplicity, we 
assume the deconfined phase is an ideal gas of free gluons and light quarks. Any J/tp 
in this medium will be subject to dissociation via collisions with gluons. (This is the 
dynamic counterpart of the plasma screening scenario, in which the color-confinement 
force is screened away in the hot dense plasma B.) The primary formation mechanism 
is just the reverse of the dissociation reaction, in which a free charm quark and 
antiquark are captured in the J/ip bound state, emitting a color octet gluon. Thus it is 
unavoidable for this model of quarkonium suppression that a corresponding mechanism 
for quarkonium production must be present. The competition between the rates of 
these reactions integrated over the lifetime of the QGP then determines the final J/ip 
population. Note that in this scenario it is impossible to separate the formation process 
from the dissociation (suppression) process. Both processes occur simultaneously, in 
contrast to the situation in which the formation only occurs at the initial times before 
the QGP is present. 

The time evolution if the J/ip population is then given by 

—£± = Xf N c P - c - XdN jh> Pg , (4) 

where r is the proper time in a comoving volume cell and pi denotes the number 
density [L~ 3 ] of species i. The reactivity A [i 3 /time] is the reaction rate {av rc \) 
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averaged over the momentum distribution of the initial participants, i.e. c and c for 
X F and J /if) and g for X D . The gluon density is determined by the equilibrium value 
in the QGP at each temperature. Exact charm conservation is enforced throughout 
the calculation. The initial volume at r = tq is allowed to undergo longitudinal 
expansion V(r) = Vqt/tq. The expansion is taken to be isentropic, VT 3 = constant, 
which then provides a generic temperature-time profile. For simplicity, we assume the 
transverse spatial distributions are uniform, and use a thermal equilibrium momentum 
distribution for both gluons and charm quarks. (This last simplification requires large 
energy loss mechanisms for the charm quarks in the deconfined medium, which is 
indicated by several recent studies ||). 

With these inputs and assumptions, the solution of Equation [| is precisely that 
anticipated in Equation |lj, with 

e( Tf ) = e J -o , (5) 

where Tf is the hadronization time determined by the initial temperature (T is a 
variable parameter) and final temperature (Tf = 150 MeV ends the deconfining phase), 
and 

0(r f ) = e(r/) x r A F [V(r) e(r)]- 1 dr. (6) 

J TO 

For our quantitative estimates, we utilize a cross section for the dissociation of 
J/tp due to collisions with gluons which is based on the operator product expansion 
y, |(J , which is utilized with detailed balance factors to calculate the primary formation 
rate for the capture of a charm and anticharm quark into the J/tp. 

4. Results 

In Figure [ll we show the time development of the J /ip (solid line) along with the 
separate formation and dissociation rates (dotted lines, arbitrary units). 

This calculation maintained exact charm conservation, so that the solutions 
followed evolution of both bound and free charm quarks. One sees the expected 
decrease of the formation rate due to the volume expansion, and the decrease of the 
gluon dissociaton rate due to the decrease in gluon density with temperature. 

Some typical calculated values of the J/ip final population are shown in Figure 0. 
The parameter values for thermalization time To = 0.5 fm, initial volume Vq = ttR 2 to 
with R = 6 fm, and a range of initial temperature 300 MeV < To < 500 MeV, 
are all compatible with expectations for a central collision at RHIC. The quadratic 
fits of Equation [l| are superimposed, verifying our expectations that the decrease 
in initial unbound charm is a small effect. (These fits also contain a small linear 
term for the cases in which N\ is nonzero, which accounts for the increase of the 
unbound charm population when dissociation occurs.) The fitted e values decrease 
quite rapidly with increasing T Q as expected, and give reasonable upper limits for 
the suppression factor of directly-produced J/ip in central collisions at RHIC due to 
gluon dissociation in a deconfined phase. The corresponding (3 values are relatively 
insensitive to To, remaining in the range 2.0 — 2.6 x 10~ 3 . These fitted parameters 
must be supplemented by values of x and No to determine the "suppression" factor 
from Equation || for the new mechanism. We use No = 10 from a pQCD estimate 
M. An order of magnitude estimate of 10 -2 for x, from fitted values of a color 
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Figure 1. Time dependence of J/ip formation including new mechanism. 
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Figure 2. Calculated J/ip formation in deconfined matter at several initial 
temperatures at RHIC, as a function of initial charm production. 



evaporation model |§|], is reduced by the suppression due to interactions with target 
and beam nucleons. For central collisions we use 0.6 for this factor, which results 
from the extrapolation of the observed nuclear effects for p-A and smaller A-B central 
interactions. 

With these parameters fixed, we predict from Equation || an enhancement factor 
for J/tp production of 3.6 < Sj/^ < 5.4, for initial temperatures between 300 and 500 
MeV. The suppression of initially- produced J/ip alone ranges from factors of 10 to 100, 
so that the enhancement prediction involves a huge increase (factors of approximately 
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one to two orders of magnitude) in the final population of J/tp at RHIC. 



5. Centrality Dependence 

We also predict how this new effect will vary with the centrality of the collision, 
which has been a key feature of deconfinement signatures analyzed at CERN SPS 
energies 0. The e and (3 parameters are recalculated, using appropriate variation 
of initial conditions with impact parameter b. From nuclear geometry and the total 
non-diffractive nucleon-nucleon cross section at RHIC energies, one can estimate the 
total number of participant nuclcons Np(b) and the corresponding density per unit 
transverse area np(b,s) [Q. The former quantity has been shown to be directly 
proportional to the total transverse energy produced in a heavy ion collision [ |ll| . The 
latter quantity is used, along with the Bjorken-model estimate of initial energy density 
[|12| , to provide an estimate of how the initial temperature of the deconfined region 
varies with impact parameter. We also use the ratio of these quantities to define 
an initial transverse area within which deconfinement is possible, thus completing 
the initial conditions needed to calculate the J/tp production and suppression. The 
average initial charm number 7Vn varies with impact parameter in proportion to the 
nuclear overlap integral Taa(\))- The impact-parameter dependence of the fraction x is 
determined by the average path length encountered by initial J/tp as they pass through 
the remaining nucleons, L(b) |13fl . All of these b-dependent effects are normalized to 
the previous values used for calculations at b = 0. 
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Figure 3. Ratio of final J/tp to initial charm as a function of centrality, 
due to nuclear absorption only (solid line), after final state suppression by 
a QGP (dashed lines), and with inclusion of the new formation mechanism 
in the deconfined medium (solid symbols). 

It is revealing to express these results in terms of the ratio of final J/tp to 
initially-produced charm pairs. In Figure 0, the solid symbols are the full results 
predicted with the inclusion of our new production mechanism at RHIC. The centrality 
dependence is represented by the total participant number iVp(b). For comparison we 
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also show predictions without the new mechanism, when only dissociation by gluons 
is included (Xp = 0). It is evident not only that the new mechanism dominates the 
J/ip production in the deconfined medium at all impact parameters, but also that 
an increase rather than a decrease is predicted for central collisions. These features 
should be distinguishable in the upcoming RHIC experiments. 

6. Model Dependence 

In our model of a deconfined region, we have used the vacuum values for masses and 
binding energy of J ftp, and assumed that the effects of deconfinement are completely 
included by the dissociation via gluon collisions. For a complementary viewpoint, we 
have also employed a deconfinement model in which the J/ip is completely dissociated 
when temperatures exceed some critical screening value T s . Below that temperature, 
the new formation mechanism will still be able to operate, and we use the same cross 
sections and kinematics. We find that for T s — 280 MeV, the final J/ip population 
is approximately unchanged, while decreasing T s to 180 MeV could reduce the J/ip 
production by factors of 2 or 3. These results are shown in Figure 0. 
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Figure 4. Comparison of Gluon Dissociation and Screening Scenarios for 
J/ip formation including new mechanism (see text for details). 



We have also checked the sensitivity of these results to several other assumptions 
and parameters. Among these are: (a) Change in initial charm production due to gluon 
shadowing; (b) Alternative cross sections with different magnitudes and threshold 
behaviors; (c) Transverse expansion of the QGP; (d) Non-chemical equilibrium for 
gluons; (e)Non-thermal momentum distributions for charm quarks. The effects of 
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varying these assumptions produce both positive and negative changes in the final 
J/ijj populations. The largest effect could be a decrease by a factor of 2 or 3 if one 
uses the initial pQCD momentum distributions for the charm quarks. Taken together, 
however, it is unlikely that a conspiracy of these effects would qualitatively change 
the predicted enhancement effects of this deconfinement scenario. 

7. Summary 

In summary, we predict that at RHIC energies the J /if) production rate will provide 
a more interesting signal for deconfinement than has been previously realized. 
Consideration of multiple heavy quark production made possible by higher collision 
energy effectively adds another dimension to the parameter space within which one 
searches for patterns of quarkonium behavior in a deconfined medium. It will be 
possible to experimentally "tune" the number of initial heavy quark pairs by sweeping 
through either centrality or energy. One can then search for a J/ip production behavior 
which is predicted to be nonlinear in total charm. In our simplified kinetic model of 
J/tp formation in a free gas of quarks and gluons, the new production mechanism 
predicts an enhancement rather than a suppression. These features should provide a 
signal at RHIC which will be difficult to imitate with conventional hadronic processes. 
The extension of this scenario to LHC energies will involve hundreds of initially- 
produced charm quark pairs and multiple bottom pairs. We expect the effects of this 
new production mechanism to be striking. 

Acknowledgments 

Acknowledgment: This work was supported by a grant from the U.S. Department of 
Energy, DE-FG03-95ER40937. 

References 

[1] Matsui T and Satz H 1986 Phys. Lett. B178 416 

[2] Schroedter M, Thews R L and Rafelski J 2000 Phys. Rev. C62 024905 

[3] Thews R, Schroedter M and Rafelski J 2000 Preprint fricp-ph/0007323 

[4] Kharzeev D and Satz H 1994 Phys. Lett. B334 155 



[5] For a review, see Baier R, Schiff D and Zakharov B 2000 Preprint frep-ph/000219£ 

[6] Peskin M E 1979 Nucl. Phys. B156, 365; Bhanot G and Peskin M E 1979 Nucl. Phys. B156 

391 
[7] McGaughey P L, Quack E, Ruuskanen P V, Vogt R and Wang X -N 1995 Int. J. Mod. Phys. 

A10 2999 
[8] Gavai R, Kharzeev D, Satz H, Schuler G, Sridhar K and Vogt R 1995 Int. J. Mod. Phys. A10 

3043 
[9] NA50 Collaboration, Abreu M C et al. 2000 Phys. Lett. B477 28 
[10] Bialas A, Bleszynski M and Czyz W 1976 Nucl. Phys. Bill 461 
[11] Margetis S et al. 1995 Nucl Phys. A590 355c 
[12] Bjorken J D 1983 Phys. Rev. D27 140 
[13] Gerschel C and Hiifner J 1992 Z. Phys. C47 171 



